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SECTION  I 


INTRODUCTION 

The  broad  objective  of  this  effort  is  to  guide  the  experi- 
mental invesitgation  of  electromagnetic  pulse  (EMP)  simulation 
by  portable  simulators.  We  perform  two  distinctly  different 
types  of  analysis  directed  toward  this  objective.  First,  we 
develop  a theory  that  results  in  the  definition  of  technical 
objectives  for  both  experiments  and  calculations.  Finally, 
we  perform  calculations  to  determine  whether  a very  idealized 
experiment  could  possibly  achieve  the  required  objectives. 

The  analysis  resulting  in  the  technical  objectives  consists 
of  developing  the  form  of  a transfer  operator  equation  in 
sufficient  detail  to  identify  the  significance  of  all  terms. 
Specifically,  attention  is  directed  toward  clearly  identifying 
the  physical  quantities  related  by  the  transfer  operator  as 
well  as  the  physical  quantities  on  which  the  transfer  operator 
depends.  To  facilitate  the  discussion  of  the  physical  quan- 
tities it  is  necessary  to  discuss  the  type  of  system  we  wish 
to  excite  with  the  portable  simulators.  The  class  of  systems 
for  which  this  study  is  applicable  are  those  systems  that  are, 
in  effect, imperfectly  sealed  metallic  enclosures.  Important 
systems  that  belong  to  this  class  are  aircraft,  missiles,  ships, 
and  tanks.  The  breaks  in  these  enclosures  are  referred  to  as 
apertures  and  they  might  correspond  to  windows,  hatches,  or 
portions  of  deliberate  antennas  that  are  intended  to  allow 
energy  to  flow  into  the  system. 

The  operator  equation  relates  electrical  quantities 
excited  within  the  actual  enclosure  (system)  to  the  current 
density  induced  on  metallic  seals  placed  over  all  of  the 
apertures  of  the  imperfectly  sealed  enclosure.  The  existence 
of  this  equation  would  seem  to  imply  that  if  a configuration 
of  portable  sources  excited  the  same  current  density  on  the 
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seals  as  did  an  EMP,  then  internal  electrical  quantities  with- 
in the  enclosure  of  the  unsealed  system  would  be  identically 
excited  by  either  the  EMP  or  the  portable  source  configuration. 
This  would  be  the  case  under  the  following  conditions;  the 
portable  sources  must  be  electromagnetically  rigid,  i.e., 
unaffected  by  the  presence  of  any  scatterer,  and  the  external 
environment  of  the  system  must  be  the  same  for  the  portable 
source  configuration  as  for  the  EMP.  For  example,  an  aircraft 
having  the  appropriate  seals  correctly  excited  by  rigid  port- 
able sources  when  parked  on  the  ground,  can  only  be  viewed 
as  having  been  excited  by  the  corresponding  EMP  when  it  is 
still  resting  on  the  ground  and  in  particular  is  not  in  free 
flight. 

Even  with  these  limitations,  we  see  that  it  is  possible 
to  assist  the  alternate  simulation  program  by  performing  only 
external  interaction  measurements  or  calculations.  The 
initial  source  configuration  can  be  determined  by  employing 
only  external  interaction  considerations.  We  emphasize  that 
we  expect  the  focus  to  be  on  external  interaction  only  in  the 
initial  program  stages  because  we  anticipate  that  the  local 
sources  will  not  be  capable  of  exciting  exactly  the  same 
external  interaction  quantities  on  the  metallic  seals  as 
would  an  EMP.  In  order  to  assess  these  effects  as  well  as 
non-rigidity  degradation,  we  expect  that  internal  electrical 
quantities  will  have  to  be  measured  for  excitation  by  the 
portable  source  configuration  as  well  as  for  excitation  by 
a more  orthodox  simulator  which  represents  the  EMP  excitation. 

The  environment  and  source  rigidity  conditions  previously 
discussed  result  from  the  dependence  of  the  transfer  operator 
on  these  factors  and  not  the  quantities  this  operator  relates. 
This  source  rigidity  requirement  causes  special  concern  in 
that  any  physically  realizable  portable  source  is  going  to 


have  structure  that  can  interact  with  the  fields  reflected 
from  the  system  under  test.  This  is  of  particular  concern 
because  it  is  presently  anticipated  that  the  configuration 
of  portable  sources  will  be  in  close  proximity  to  the  system 
under  test.  The  choice  of  calculations  to  perform,  which 
represented  an  idealized  experiment,  was  made  with  the  issue 
of  source  rigidity  being  a distinct  factor. 

The  problem  for  which  we  made  our  calculations  was  the 
excitation  of  a sphere  in  free  space  by  a plane  wave  and  by 
various  conf igurations  of  idealized  local  sources.  These 
calculations  were  performed  in  the  frequency  domain  for  a 
range  of  frequencies  starting  at  zero  and  extending  to  approx- 
imately three  times  the  first  resonant  frequency  of  the  sphere. 
Well-established  pliine  wave  solutions  exist  for  this  problem 
and  our  method  of  obtaining  our  plane  wave  solution  can  be 
verified  by  comparison  of  our  results  to  the  established  re- 
sults. This  is  necessary  because  our  method  of  obtaining  the 
plane  wave  solution  is  the  same  as  our  method  of  obtaining  the 
source  configuration  results  and  no  data  is  presently  available 
to  verify  those  calculations.  As  a general  conclusion,  our 
calculations  indicate  that  our  choice  of  local  source  config- 
uration can  approximately  excite  the  desired  external  inter- 
action current  density  at  a shorted  point  of  entry  only  if  at 
least  one  local  source  is  in  close  proximity  to  the  shorting 
surface.  This  result  increases  the  need  to  study  the  effect 
of  the  degree  of  rigidity  of  physically  realizable  sources  on 
the  alternate  simulation  problem. 


SECTION  II 


THEORETICAL  BACKGROUND 


This  investigation  concerns  the  local  excitation  of  systems 
that  are  predominantly  metallic  and  is  valid  for  those  fre- 
quencies or  times  for  which  the  metal  can  be  considered  to  be 
perfectly  conducting.  The  equations  that  form  the  basis  of 
this  investigation  of  portable  EMP  simulators  is  a set  of 
equations  that  recognizes  those  essential  features  of  classical 
aperture  coupling  analysis  that  have  relevance  to  complex 
systems.  Since  this  approach  is  based  on  aperture  coupling 
equations,  one  might  be  concerned  with  its  relevance  to  other 
types  of  penetrators,  e.g.,  deliberate  antennas.  Such  pene- 
trators  have  associated  apertures  or  else  no  energy  could 
penetrate  the  sealed  skin  of  the  system  corresponding  to  that 
penetrator . 


First  we  will  present  the  general  form  of  the  equations 
that  provide  the  basis  of  this  study  and  draw  all  of  our 
theoretical  conclusions  by  referring  to  properties  of  this 
general  form.  Next  we  will  present  a somewhat  detailed 
derivation  of  these  general  equations  for  a complex  interaction 
situation  in  order  to  give  a more  concrete  meaning  to  the 
general  properties  on  which  we  based  our  theoretical  conclusions. 
The  form  of  the  underlying  equation  is  as  follows 


= — E . 


I.(£> 


(1) 


where  the  meaning  and  significance  of  each  term  requires  consider- 
able attention.  First,  we  emphasize  that  equation  1 describes 
the  relationship  between  electrical  quantities  on  two  different 
physical  systems.  One  system  is  the  actual  system  of  interest 
and  the  other  system  is  that  original  system  modified  by 
metallic  shorting  surfaces  covering  all  apertures  (including 
those  associated  with  antennas).  For  illustrative  purposes 
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consider  the  system  depicted  in  figure  1.  One  system  is  the 

aircraft  in  its  environment  with  the  apertures  and 

unmodified  and  the  other  system  needed  to  give  equation  1 

meaning  is  the  same  aircraft  in  the  same  environment  with 

metallic  seals  covering  and  S0.  In  equation  1,  the 

notation  J was  chosen  to  denote  "magnetic  current,"  but 
— m 

it  is  simply  n (r ' ) x E (r ' ) where  r'  varies  over  all  of  the 
mathematical  surfaces  corresponding  to  the  open  apertures  in 
the  original  system,  n (r ' ) is  the  outward  normal  at  r ' , and 
(r ' ) is  the  tangential  component  of  the  electric  field 
induced  in  the  open  aperture.  The  quantity  J„,T  (r)  is  the 
"external  interaction"  current  density  induced  on  the  shorted 
system  with  r ranging  only  over  the  shorting  surfaces.  It  is 
important  to  note  that  even  though  r and  r'  refer  to  different 
physical  systems,  they  mathematically  refer  to  the  same  set 
of  points.  This  distinction  allows  a discussion  of  the  mathe- 
matical nature  of  equation  1 that  is  not  confused  by  the 
dual  physical  nature  of  the  problem.  It  remains  to  discuss 
the  meaning  of  L in  equation  1 to  proceed.  More  specifically, 
L is  a linear  operator  that  depends  on  a variety  of  quantities 
associated  with  the  system,  its  environment  and  certain 
aspects  of  its  excitation.  Just  what  these  quantities  are 
plays  an  essential  role  in  the  underlying  theory  of  portable 
EMP  simulators  and  we  will  elaborate  on  what  these  quantities 
are  when  presenting  the  details  for  the  system  depicted  in 
figure  1. 

It  is  now  necessary  to  introduce  an  additional  equation 
to  augment  the  information  contained  in  equation  1.  This 
equation  also  represents  a general  form  and  is 

(2) 

This  equation  is  a mathematical  statement  of  the  fact  that 

J (r* ) is  sufficient  to  determine  a variety  of  electrical 

— tn  — 
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Q0(e.g.,  8 can  correspond  to  a voltage,  a current,  or  a field 
p 

component)  that  are  excited  within  the  system  by  fields 

penetrating  through  the  apertures.  In  equation  2,  L„  is  a 

P 

linear  operator  that  depends  on  the  internal  structure  of  the 
system  and  the  choice  of  the  internal  electrical  quantity 
that  is  heing  determined.  Next  we  introduce  a step,  the 
legitimacy  of  which  is  currently  being  studied  using  a field 
equivalence  point  of  view.  Specifically,  it  is  assumed  that 
the  L appearing  in  equation  1 has  a unique  inverse,  L so 
that  from  equation  1 we  can  obtain 


(3) 


Combining  equations  2 and  3 we  obtain 


where 


°s 


I. 


(4) 


tb  = V'1  <5) 

and  the  superscript  a is  explicity  introduced  to  indicate 

that  depends  on  the  environment  external  to  the  system. 
p 

If  the  same  system  were  placed  in  two  different  environments, 
then  the  a designation  for  each  environment  could  change  to 
accommodate  a mathematical  representation  of  the  fact  that 

a.  a 

T6  + T62  (6) 

if  the  external  environments  for  the  same  system  are  sufficiently 
different.  Part  of  what  we  shall  mean  by  the  external  environ- 
ment is  the  physical  structure  of  the  portable  EMP  simulators 
that  are  being  investigated.  When  we  discuss  the  details 
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with  the  system  and  environment  depicted  in  figure  1,  we  shall 
emphasize  this  source  structure  dependence  and  make  a crucial 
distinction  between  rigid  and  nonrigid  sources. 

It  is  possible  to  present  all  of  the  portable  simulator 
theory  on  equation  4;  however,  that  equation  will  be  modified 
to  conform  to  the  prevalent  notion  that  both  the  external 
interaction  current  density  J and  the  external  interaction 

Cj  • 1 • 

charge  density  ap  T are  required  for  the  ultimate  determine- 

tion  of  the  internal  quantities  Q„.  For  non-zero  frequency, 

p 

it  follows  from  V *J  = iwo  that  J suffices  to 

S lii  • X • lii  • X • ill  • X • 

determine  o , so  the  requirement  that  o_  be  separately 
Xi  • X • E • X • 

determined  must  be  superfluous.  There  are  a number  of  pos- 
sibilities why  it  might  be  convenient  to  separately  view 
o„  T as  a desired  input  and  viewing  it  as  such  leads  to  the 
following  decomposition  of  equation  4 


QS  TJB-E.I.  + To  8°E . I . 


(7) 


as  the  basic  equation. 

At  this  point  we  could  present  the  underlying  theory  of 
portable  EMP  simulators  by  referring  to  either  equation  4 
or  equation  7 if  we  did  not  have  to  deal  with  the  real 
physical  structure  of  the  portable  sources. 

The  means  whereby  this  aspect  enters  the  consideration  is 
rather  complex  and  is  treated  by  giving  a more  explicit  meaning 
to  these  equations.  Specifically,  this  will  be  accomplished 
by  deriving  more  explicit  representations  for  equation  1 
and  equation  2 for  the  situation  depicted  in  figure  1.  First 
we  introduce  the  following  definitions: 

S : the  surface  of  the  metallic  enclosure  (aircraft) 
m 

augmented  by  the  mathematical  surfaces  and 

V^:  is  the  volume  of  a lossy  medium  in  the  proximity  of 

the  enclosure  (earth,  water) 
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the  surface  bounding  VL 

the  volume  of  an  object  in  the  proximity  of  the 
enclosure  (i.e.,  an  aircraft  carrier) 
the  surface  hounding  V 

P 

the  volume  of  a subsystem  contained  within  the 
enclosure 

the  surface  bounding  V 

the  volume  exterior  to  S bounded  by  S , S , St  , 

m m p L 

S^,  and  the  hemisphere  at  infinity 

the  volume  interior  to  S bounded  by  Sm  and 

m -‘ms 

the  volume  of  a rigid  source  of  an  electromagnetic 

wave,  J,  and  it  is  contained  in 

the  volume  of  the  portable  radiator 

the  surface  of  the  portable  radiator 

the  portion  of  S over  which  the  surface  tangential 

electric  field  is  rigidly  specified 


The  essential  equation  that  this  approach  is  based  on  is  the 
dyadic  identity 


A (r ' ) • | V ' xV ' xD (r ' ;r)j  -|v'xV'xA(r'  ) J • D (r 1 ; r ) 


dV' 


A (r ' ) 


n (r' ) x VxD  (r’ 


+ |(v  1 xA  (r 1 ))  xn  (r ' ) ] • D (r ' ; r 

(8) 


where  A(r ' ) and  g (r ' ;r)  are,  at  this  point,  a general  vector 
and  a general  dyadic  that  must  satisfy  certain  behavior 
requirements  (e.g.,  differentiability)  but  not  necessarily 
any  equations.  In  equation  8,  S is  the  surface  bounding  V 
and  n (r ' ) is  the  outward  normal  to  V.  Next,  the  volume, 
bounding  surface,  A(r' ),  and  g(r';r)  are  specialized.  V is 
chosen,  in  turn,  as  VQ  and  and  A(r' ) is  chosen  as  (r ’ ) 
and  ^ (r ' ) . We  also  choose  D (r ’ ;r)  as  appropriate  Green's 
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- 


dyadics  ( (r' ; r) , gj  (£';£) j that  satisfy  the  vector  wave 
equation 


(v'xVx-kgj  G 


(r'  , 
a -a 


= I6(r'-r  ! 
= —a  —a 


a*0 . I 


r (V 


(9) 


and  subsequently  the  a subscript  of  r and  r'will  automatically 
be  implied  by  the  subscript  on  when  it  is  not  explicitly 
indicated.  Boundary  conditions  to  be  satisfied  are 

n(r'  )x(v,xGJ  (r'  ,r))  = 0 r'<Sm  (10) 


n (r* ) x(v' xgQ  (r* ,r))  = 0 r’<SrUSpUSm  (11) 

n(r')x(7'xg  (r',r))»  n (r ' ) x(  V ' xgL  (r ' , r ))  r'fSL  (12a) 


n (r ' ) xe0g0 (r ' ,r)  = n (r ' ) xegL (r ' , r ) r'f  SL  (12b) 


The  equation  satisfied  by  G^(r',r)  is 

xV' x-u)2uQc.  jgL(r' ,r)  = 0 r'fV  j ,r£VQ 


(13) 


The  equations  satisfied  by  the  H((r')  are 

V ' XV 1 x-k?  In  (r  ■ > » j 0 <'”1 

0,-“  " | VxJ(r')  a=0 

It  also  follows  from  Maxwell's  equations 


(14) 


VxH^r'l  = -lue^fr')  r'!S<  (15) 

Substituting  equations  9,  14,  and  15  into  8 for  V=Vq  or 
and  using  the  property  of  the  <$  function,  we  obtain 
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H0(I0)  = Kr 


0>  ^ J |nq'-0'-0'=0 

q=m, p, r , L S | 

+I  |v 


dS1 


'lo^o'So  ds’ 


«i<Si>  - XI  / jv-i'-i'£i  j ds' 


q=m,s  S 


(16) 


(17) 


where  afI  is  the  unit  outward  normal  to  the  sphere  at  infinity 

S.  , 


/ j n , E , H , Q J 


dS' 


1 

LI  11 

A 

1 

j H(r’  ) • 

n (r 1 ) x V' xQ (r'  ; r)  j 

+ iwe 

n (r  ' ) xE  (r  ' ) 

•Q (r ' ; r)  j 

where  c is  the  appropriate  dielectric  permittivity  and 


V'xJ  (r' ) *g0 (r' ;rQ)  dv' 


(18) 


(19) 


Using  equations  10  and  11  as  well  as  the  fact  that 


n (r  * ) xE  (r  ' ) = 0 r'f(S  -S  -S.)US  U (S  -S  ) 
— — — — m 1 2 p r g 


we  find  that 


f 

1 

n , E.  , 11  , G - 1 

' S 

P 

-0-0=0 j 

(20) 


(21) 


dS' 
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Is 


g iwfQ  |n  (r'  )XE  (r ' )J  *GQ  (r ' /r^)  dS' 
9 


5 do1 


(22) 


i ioe 


“(iK- 


II  «■ 


' ) XE  (r* 
—a  — 


■G  (r’;r  ) dS' 
= a — —a 


) XE^  (r ' ) 


dS- 


a = 0,1 


(23) 


and  because  EQ,  H , and  gQ  satisfy  the  radiation  condition 


(24) 


The  remaining  quantities  to  evaluate  in  equations  16  and  17  are 

the  surface  integrals  over  ST  and  S . Substituting  the  equa- 

lj  s 

tions  appropriate  for  the  lossy  half  space,  that  is 


(V*  xV'  x-w‘ 


MqOHlU') 


= 0 r'(V, 


(25) 


and 


V'xHL(r')  =-iweEL(r')  £,fVL 


as  well  as  equation  13 


into  8 we  obtain 


(26) 
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f 

Ht  (r*  ) • 

n (r ' ) x 

V*  XQ  (r ' ; r ) j 

+ iuie 

n (r ' )XE  (r’  ) 

G.  (r'  ; r) 

/SL 

— L — 

~L  ~ ~ IJ 

— Li 

dS1 


/• 

1 r 

( 

oo 

HL(r' ) • 

n (r ' ) x j V ' xGl (r ' ;r) J 

+ iwe 

n (r1 )xEr  (r ' ) 

Li 

GL( r' ;r) 

dS1 


(27) 


The  second  integral  in  equation  27  is  zero  due  to  the  losses 
in  (or  the  radiation  condition  if  VL  is  lossless).  Using 
the  fact  that  the  tangential  components  of  E and  H are  contin- 
uous across  S as  well  the  boundary  conditions  in  equation  12a 
and  12b  we  see  that  the  integral  over  S in  equation  16  is  equal 

L« 

to  the  integral  over  S in  equation  27  which  in  turn  we  have 
just  shown  to  equal  zero.  The  integral  over  Sg  will  also 
equal  zero  and  the  manner  in  which  this  can  be  seen  depends 
on  the  physical  properties  of  the  subsystem  occupying  V . If 
it  were  totally  metallic,  the  boundary  conditions  on  E and 
G would  make  the  surface  integral  vanish  in  the  same  manner 
they  did  for  the  integral  over  S^.  If  it  were  a homogeneous 
dielectric,  then  the  boundary  conditions  would  cause  the 
surface  integral  over  Sg  in  the  same  manner  the  surface 
integral  over  S was  caused  to  vanish.  If  it  were  some 
hybrid  of  dielectric  and  metal,  a combination  of  the  arguments 
would  be  used  to  cause  the  surface  integral  to  vanish. 

We  can  now  write  equations  16  and  17  as 


!W  ’ - WE1) 


(28) 


and 


-I (-I J = 


(29) 
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where 


*<£ o>  " 


+ i^0) 


(30) 


with  I(£q)  and  S(rQ)  defined  by  equations  19  and  22  and 
operators  are  defined  by 


K J (r‘)  = iu>en(  / J (r')«G  (r',r 
at— m — 0^  J g — m — =cx  — ct 


Is. 


) dS' 


+ I J (r* ) *G  (r' ,r  ) dS' 
* — m =a  — — a 


the 


a = 0,1  (31) 

and  we  have  made  use  of  the  fact  that  the  tangential  components 
of  the  electric  field  are  continuous  through  the  apertures  so 
that 


-n0(r'  )xE0(r'  ) = nx  (r  • ) xEj  (r ' ) = Jjr’  ) r'eS1US2  (32) 

Now  we  focus  our  attention  on  FCr^)  appearing  in  equation  28. 
The  meaning  of  this  quantity  is  an  extremely  important  aspect 
of  the  theory  behind  portable  EMP  simulators.  It  would  be 
a very  difficult  task  to  evaluate  equations  19,  22,  and  30 
in  order  to  determine  the  full  significance  of  F(r^).  Instead, 
we  will  simply  utilize  certain  key  features  of  those  equations 
as  well  as  equations  28  and  31  to  determine  what  F(£q)  must 
be  if  all  the  required  equations  were  evaluated.  First,  we 
note  according  to  equation  19  that  £(rQ)  is  excited  by  the 
rigid  (interaction  independent)  source  J(r' ) and  that  accord- 
ing to  equation  22,  Sfr^)  is  excited  by  the  rigidly  specified 
n(r')xE(r')  for  r'fSg.  Next,  we  note  that  according  to  these 
equations,  both  Mr^)  and  S(£g)  are  insensitive  to  the  size 


of  the  apertures  and  S2  and  in  fact  they  are  insensitive 
to  whether  or  not  these  apertures  are  even  present.  Using 
these  observations  in  conjunction  with  equation  31  as  the 
aperture  size  becomes  zero  and  using  the  result  in  equation 
28,  we  see  that  F(£q)  equals  (r)  for  the  special  case 
where  all  apertures  are  sealed  (short  circuited) . Mathema- 
tically, we  express  this  evaluation  of  F(£q)  as 


«o*C*  (^0> 


(33) 


where  the  superscript  is  introduced  to  indicate  "short 
circuit."  We  note  that  £(£q)  is  the  short  circuit  magnetic 
field  at  some  point  £q  with  apertures  sealed,  but  all  other 
aspects  of  the  external  environment  including  the  proximity 
and  structure  of  the  radiator,  Sf,  unchanged. 

Substituting  equation  33  into  equation  28  we  obtain 


VE0>  * 


(34) 


A A ^ 

Next,  we  define  n (£)  = n^.  (£)  = -nQ  (£)  for  rfS^US^  use 
the  fact  that 


n (£)  XHq  (£)  = n(£)xH];(£) 
and  employ  equations  28  and  34  to  obtain 


(35) 


n(r)x  (KQ+KI)Jn(r' ) =JE>I>(£)  (36) 


where  we  have  used  the  definition 

n(£)XHg-C- (£)  = (£)  (37) 


lim 

-I**- 
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r 


t 


1 


i 


and  we  have  the  desired  result,  in  that  equation  36  is  the 
more  detailed  representation  of  equation  1. 

Before  we  can  present  our  theoretical  conclusions , we 

must  present  our  more  detailed  representation  of  equation  2. 

We  have,  in  fact,  already  a representation  of  equation  2 for 

the  case  where  the  desired  internal  electrical  quantity  is  the 

magnetic  field.  For  that  use  we  might  choose  the  symbol  6 

as  H so  that  Q =H  and  L =L  =K_.  Another  example  where  the 
H —■*  p n X 

structure  of  Lg  changes  depending  on  the  choice  of  Qg  is 
readily  demonstrated  by  considering  the  case  where  the  desired 
internal  electrical  quantity  is  the  electric  field  E and  we 
denote  8 as  E so  that  QE=E.  For  this  case  equation  2 becomes 

0E  - Lgjr.Or)  (38) 


where 


lL>£ 


-VxK, 


(39) 


Finally,  we  will  discuss  the  more  important  case  where  the 
desired  internal  electrical  quantity  is  a current.  For  this 
discussion  consider  that  part  of  the  internal  subsystem 
occupying  volume  V in  figure  1 contains  a wire  and  we  choose 

S 

a local  cylindrical  coordinate  system  having  its  axis  along 
the  wire  and  having  the  local  azimuthal  vector  denoted 

at  the  point  on  the  wire  where  we  wish  to  determine 
the  current.  The  argument  of  this  unit  vector,  l' , denotes 
the  circumferential  position  on  the  wire.  With  these  defini- 
tions, the  current  on  the  wire  is 


I 


(40) 
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We  see  that  from  equations  29  and  40  that 


Q = L J (r' 
c c-m  — 


where  we  have  denoted  I = Q and 


L J ( 
c-m 


r'  ) s J dr 


Vr|-KIV£') 


We  have  now  presented  equations  1 and  2 in  sufficient 
detail  to  draw  our  desired  conclusions.  We  will  base  our 
conclusions  on  equation  4 which  contains  exactly  the  physics 
as  do  equations  1 and  2.  The  specific  points  we  wish  to 
make  are  i)  the  external  interaction  current  density,  J£  , 
can  be  excited  by  either  a rigid  source,  a non-rigid  source, 
or  a combination  of  the  two  types  ii)  the  transfer  operator, 

Tg,  depends  on  the  external  environment  to  the  system  iii) 

a q 

T„  depends  on  the  internal  environment  iv)  Tg  depends  on  the 

internal  electrical  quantity,  Qe,  being  determined  v)  T„ 

p p 

depends  on  the  rigidity  of  the  source.  Equations  that  speci- 
fically illustrate  each  of  these  points  are  identified  with  the 
numbered  points  as  follows:  i)  equations  19,  22,  30,  33,  and 
37  ii)  equations  11  and  12  iii)  equation  10  as  well  as  the  argu- 
ment that  eliminated  the  integral  over  Sg  iv)  equations  5,  39, 
and  42  v)  equation  11. 

The  remaining  portion  of  this  report  will  be  devoted  to  the 
calculation  that  represents  the  idealized  experiment. 


SECTION  III 


MAGNETIC  FIELD  INTEGRAL  EQUATION  FOR  A SPHERE 

_ A A 

II  we  impose  an  orthonormal  coordinate  system  s,t  on  a 
closed  surface  possessing  continuous  curvature  such  that 
sxt  = n, the  outward  normal  to  the  body,  we  can  write  the 
Magnetic  Field  Integral  Equation  (MFIE)  as  the  following 
system  of  coupled  scalar  integral  equations 


f Js{I)  = 


-t  (r) Mllnc  (r)  + 


A(r,r' ) Jg  (r'  ) 


+ 13  (£,£'  ) Jt  ( £ 1 ) J ds' 


i nr 

1 J (£)  = S (£)  Mi  (£)  + 

2 L 


(c(£,£- 


JJs(£,) 


+D (£,£' ) J (£' ) | dS' 


where 


A(£,£’ ) 

= -Q  ( 1 £“£’  | ) | 

t (£)  ’ (£“£'  ) XS  ( £ ' ) 

A /V 

B ( r , r ' ) 

= -Q( | r-r’ | ) 

t (r ) • (r-r')xt(r') 

| 

N A 

C(£,£’  ) 

= Q ( | £-£'  | ) j: 

3 (£) * (£-£' ) XS (r ' ) 

D(£,£'  ) 

= Q ( | £-£'  | ) 1 S (£)•(£-£')  Xt  (£')  ] 

ikR 

Q(R)  = (ikR— 1)— T 

4i'  R 


and  J (r) , J,  (r)  are  defined  through 
s — t — 

J(£)  = Jg  (£)  s (£)  + J (£)  t (£) 


(43a) 


(43b) 


(44a) 
(44b) 
(44c) 
( 4 4d ) 

(45) 


(4b) 


A 

For  a sphere  of  radius  "a"  centered  at  the  origin,  n=(r/a) 
for  all  points  on  the  surface.  This  permits  us  to  greatly 
simplify  the  form  of  equations  44a  - 44d  even  before  specifying 
our  actual  choice  for  s and  t.  Formal  manipulation  of  the 
triple  products  in  these  equations  yield 


A (r , r ' ) = -aQ ( r- r ' ) 


s (r ) • s (r ' ) -t (r) • t (r ' ) 


B(r,r')  = -aQ(r-r')  s (r) • L (r ' )+t  (r) *s  (r  1 ) 


C ( r , r ' ) = aQ  (r-r ' ) -t (r ) • s (r ' ) -s (r ) • t (r ' ) 


(47a) 

(47b) 

(47c) 


and 


D (r , r ' ) = aQ  (r-r ' ) 


-t(r)  • t ( r 1 )+s(r)  • s ( r ' ) 


( 4 7d ) 


thus  showing  that  D(r,r')  = -A(r, r’)  and  C(r,r')  = B(r,r'). 

We  also  note  that  A(r,r')  = A(r',r)  and  B(r,r')  = B(r’,r). 
This  latter  symmetry  property  is  of  considerable  importance 
for  analytic  treatments  of  the  MFIE  on  a sphere,  but  will  be 
lost  in  the  numerical  scheme  for  solving  the  equations. 

A numerical  implementation  of  equations  43a  and  43b  even 
with  the  simplifications  of  equations  47a  - 47d  requires  that 
the  sphere  be  imbedded  in  some  coordinate  system.  We  use 
a spherical  coordinate  system,  i.e.,  an  arbitrary  position 
on  the  surface  of  the  sphere  has  cartesian  coordinates 


r(9?4>)  = a (cos<£sin0  , sin<£sinO  , cosO  ) 


(48) 


We  may  define 
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t(0,<j>)  = -eQ  = — g-g-  r(9 , 4> ) = C-cosij)cos0  , -sin4>cos0  , sin0  ) (49a) 

and 

S(6,0)  E " as£nQ  1^-  £(d  ,<p)  = (-sirup  , cos<p , 0 ) (49b) 

obtaining 

n(9,^)  = s (9  , <|>)x  t (9  , $ ) 5 ^-r(9,<J>)  = (cos4>sin0,  sin<j>sin0  ,cos6 ) 

a 

(49c) 

as  it  should. 

Inserting  equations  48,  49a  and  49b  into  equations 

47a  ~ 47d  and  recalling  that  the  element  of  area  on  the  surface 
2 

of  a sphere  is  a sin0d0d$  completes  the  specification  of  the 
MFIE  in  a spherical  coordinate  system. 

Our  procedure  for  solving  the  coupled  scalar  equations  of 
the  MFIE  is  to  partition  the  sphere  into  zones  S1  by  an 
algorithm  which  has  the  maximum  separation  of  any  two  points 
of  any  zone  tend  to  zero  as  the  number  of  zones  tends  to 
infinity.  We  then  approximate  both  Jg  and  J by  piecewise 
constant  functions  whose  discontinuities  occur  at  the  zone 
boundaries.  If  we  pick  a representative  point  from  each  zone 
and  restrict  r to  this  set  of  points,  we  obtain,  as  a matrix 
approximation  to  the  MFIE, 

dS' 


dS'  (50a) 
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I = s(r^ 


) dS' 


i ) • Hinc  (r  i ) +2  Js(r.)f  C(  £i,r' 

3 •'Sj 

♦EwI  D (r i , r ’ ) dS1 


:50b) 


This  method  of  solution  can  be  viewed  as  either  a method  of 
moments  solution  or  as  a product  integration  method. 

We  must,  however,  consider  the  nature  of  the  integrands 
in  equations  50a  and  50b.  One  can  show  that  for  an  arbitrary 
body  with  everywhere  continuous,  non-zero,  local  curvature 
A,  D,  C,  and  D are  singular  but  behave  at  worst  as 
a/|r-r‘ | as  r’  approaches  r for  some  finite  u.  This  will  be 
explicitly  shown  for  the  case  of  a sphere.  Since  we  are 
dealing  with  a two-dimensional  integral,  these  integrands 
are  still  absolutely  integrable,  however,  these  singularities 
should  be  treated  analytically  in  order  to  avoid  convergence 
problems  for  numerical  integration.  Our  programs  for  scat- 
tering from  cylindrical  bodies  remove  this  singularity  before 
attempting  the  numerical  integration;  experience  indicates 
that  such  treatment  greatly  improves  the  accuracy  of  both  phase 
calculations  and  resonant  phenomena. 

For  a sphere,  the  numerical  problem  is  much  simpler.  As 
will  be  shown  by  the  following  analysis,  a symmetric  integration 
procedure  will  permit  the  singularities  to  be  ignored  for 
sufficiently  large  zones.  By  expanding  the  scalar  triple 
products  to  second  order  in  0-0'  and  we  will  show,  as  we 

mentioned  earlier,  that  the  above  mentioned  singularity  does 
exist,  but  numerical  techniques  exist  which  avoid  the  need 
to  treat  the  singularity  analytically. 

2 

We  start  by  expanding  |r-r'  | in  powers  of  (0-0')  and  ($-<}>'). 

Ir-r'l2  = I r ' I 2 + Irl2  - 2(r-r')  = 2a2(l-n-n')  (51) 
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which  by  equation  49c  yields 


|r-r|2  = 2a2  |l-sin0sin9'(cos4>cos<}> ' +sin<{>sin<J> ' ) -cos9cos9  'J 

= a2  [sin20  (<(>-$' )2  + (0-6' )2]  + o [(9-0' )2  + (4>~4>  ’ ) 2 ) (52) 

Similarly,  from  equations  49a,  49b  and  47a  we  get 
-A (r,  r ' ) 

— ^0 = sin(j>sin4>'  + cos<|>cos4>' 

- (cos<}>cos<f> ' cos0cos0  ' +sin<fcsin<}> ' cos0co"s0  '+sin0sin9'  ) 


1 

2 


(9-0'  )2-sin20  (<*»-<*> ' ) 2|+o 


(0-4> ' ) 2+  (9-e  • ) 2 


(53) 


while  equations  49a,  49b  and  47b  yield 

-B ( r , r ) f 2 2l 

— ^ = -a  (0-0  ' ) (0-0'  ) sin©  +0  [(0-0  * ) +(9-9'  ) | (54) 

The  above  analysis  has  shown  that  neither  A nor  B behaves 
any  worse  than  a/|r-r'|  for  some  finite  a,  yet,  except  at  the 
poles  (B  is  non-singular  if  0=0)  there  exists  directions  of 
approach  such  that  both  A and  B vary  as  l/|r-r'  j as  r' 
approaches  r.  In  addition,  we  have  shown  that  except  at 
0=0  and  0=tt  B is  antisymmetric  in  (6-0')  and  ( 0— 0 ' ) and  A 
is  antisymmetric  in  (0-0')  + (0—0 ' ) sin0.  Thus  if  our 
integration  scheme  is  symmetric  in  (0-0')  and  (0—0 ' ) 

+ ($-♦') sin0  the  singular  part  effectively  vanishes  for  self 
term  interactions,  i.e.,  when  i=j  for  equations  50a  and 
50b.  However,  neighboring  zone  interactions  do  not  necessarily 
have  this  antisymmetry  property.  If  wavelength  considerations 
force  the  zones  to  be  small  the  singularities  should  be  treated 
analytically. 


Our  experience  has  shown  that  the  zoning  criteria  for 
accurate  solution  of  the  MFIE  can  be  split  into  wavelength  and 
geometry  considerations.  As  a general  rule,  between  six  and 
ten  zones  per  wavelength  are  needed  to  fulfill  the  wavelength 
requirements.  For  this  special  case  we  found  that  we  could 
employ  even  fewer  zones  for  wavelength.  For  low  frequency, 
however,  geometric  considerations  dominate  the  zoning -criteria. 
The  adequacy  of  the  geometric  requirements  can  be  ascertained 
by  examining  the  results  for  magnetostatic  excitation.  Study- 
ing both  types  of  zone  requirements,  we  found  that  the  nearest 
neighbor  zones  are  far  enough  removed  to  permit  simple  inte- 
gration schemes  for  evaluating  the  integrals  of  equations  50a 
and  50b. 
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SECTION  IV 


PRESENTATION  OF  SPHERE  CALCULATIONS 

The  coordinate  system,  incident  field  description,  and 
zone  numbering  scheme  for  this  calculation  are  depicted  in 
figure  2.  The  boundaries  for  each  zone  are  determined  by 
allowing  45°  increments  in  0 and  In  figures  3 through  11 
we  present  the  current  density  induced  by  the  depicted  incident 
field  as  well  as  by  selected  local  excitation.  What  is  meant 
by  the  local  excitation  is  that  a numbered  patch  is  either 
considered  to  be  illuminated  by  the  depicted  incident  field 
or  is  considered  to  receive  no  incident  illumination.  A 
discussion  of  the  relevance  of  this  type  of  local  illumination 
will  be  deferred  to  the  next  section. 

The  labeling  of  the  tangential  components  of  the  induced 
current  density  is  as  follows 

Js  = J<t>  <55) 

Jt  = -J9  (56) 

and  the  quantities  plotted  are  the  magnitudes  of  these  compo- 
nents of  current  density  normalized  to  the  magnitude  of  the 
incident  magnetic  field,  HQ.  The  code  verification  data 
presented  in  these  figures  comes  from  two  sources.  For  ka=0, 
the  magnitude  of  the  magnetostatic  solution  given  by 

Jg  - -(3/2)  HqcosiJ)COs8  (57) 


and 


Jt  - -(3/2)HQsin$ 


(58) 
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is  used  to  obtain  the  code  verification  data.  For  k =1.1,  1.7, 
2.3,  2.9  we  use  the  data  presented  in  figure  66a  of  reference  1. 
Specifically,  we  relate  their  data.  Kg  and  to  the  code 

verification  data  using  the  relations 


and 


jg(9,0) 

Ho 

Jfc (9,0) 

Ho 

K (y,  ) cos0 
z 1 


Kg  (6^)  sin<}> 


(59) 


(60) 


as  well  as  making  the  identification  4'^=6^=9.  The  values  of 

9 and  0 which  are  chosen  for  the  evaluation  of  equations 

57,  58,  59,  and  60  for  the  code  verification  data  correspond  to 


the  angular  centers  of  the  patches.  Finally,  we  note  that 
we  need  only  present  our  incident  field  results  for  zones 
1 through  8 because’  those  results  can  be  translated  to  the 
remaining  range  of  0 values  through  the  relations 


and 


Js(e,*p) 

H0COS*p 


Jt  (9,0) 

Jt(9,0p) 

Ho 

HoSin0p 

jg  0,0) 


H, 


COS  0 


sin0 


(61) 


(62) 


where  0^  corresponds  to 


for  zones  1 through  8. 


a value  of  0 in  the  data  presented 


1.  King,  R.W.P.  and  T.T.  Wu,  The  Scattering  and  Diffraction 
of  Waves,  Harvard  University  Press,  Cambridge, 
Massachusetts,  1959. 
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Figure  3.  Normalized  Current  Densities  on  Patch  1 
Selective  Patch  Excitation  Is  Achieved  by 
Exciting  all  Patches  but  No.  1. 


Figure  9,  Normalized  Current  Densities  on  Patch  7 
Selective  Patch  Excitation  is  Obtained  by 
Exciting  all  Patches  but 
Nos.  2,  3,  14,  15,  18,  19,  30,  31. 


SECTION  V 


mm 


INTERPRETATION  OF  SPHERE  CALCULATIONS 

The  basis  for  our  choice  of  the  illumination  scheme  that 
was  used  to  obtain  our  data  is  as  follows:  i)  the  objective 
of  each  portable  simulator  was  determined  only  by  the  incident 
field  ii)  it  was  easy  to  numerically  implement  iii)  it  bore  a 
relation  to  an  identifiable  class  of  real  sources  iv)  it  had 
to  succeed  as  more  sources  where  included.  The  choice  of  where 
to  place  the  sources  is  related  to  the  source  rigidity  issue. 
This  is  readily  seen  by  interpreting  the  results  presented  in 
figures  3 and  4.  In  each  of  these  figures,  31  of  the  32 
patches  were  illuminated  in  exactly  the  same  manner  that  they 
would  be  by  the  incident  plane  wave.  The  only  patch  that  wasn't 
excited  is  the  patch  on  which  we  present  the  data  and  we  see 
that  the  induced  current  density  is  a very  poor  approximation 
to  the  desired  current  density  which  was  induced  by  the  inci- 
dent plane  wave.  This  implies  that  if  the  non-il luminated 
patch  corresponds  to  the  shorted  POE  location,  we  can  obtain 
good  excitation  of  that  POE  only  by  having  a source,  of  the 
type  considered  in  this  report,  in  close  proximity. 

This  result  enhances  the  importance  of  source  rigidity 
effects.  This  ic  the  case  because  a qualitative  examination 
of  the  equations  that  raised  the  issue  of  source  rigidity  indi- 
cates that  the  nonrigidity  effect  becomes  increasingly  impor- 
tant as  the  source  location  approaches  the  POE.  Determining 
the  quantitative  effect  of  source  rigidity  appears  to  be  an 
experimental  problem.  Figures  7 and  8 show  that  patches  adja- 
cent to  the  nonexcited  patch  can  be  excited  in  the  desired 
manner  for  the  described  31  out  of  the  32  patch  illumination. 
This  result  again,  is  only  meaningful  if  source  rigidity  is 
not  found  to  be  a limiting  consideration. 
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The  remaining  data  appearing  in  figures  5,  6,  9,  and  10 
correspond  to  an  illumination  scheme  in  which  only  24  of  the  32 
patches  are  illuminated.  The  basis  for  choosing  not  to  illu- 
minate the  eight  patches  is  that  they  correspond  to  the  small- 
~ i nc 

est  values  of  nxii  . For  this  more  sparce  illumination  scheme 
we  again  see  that  we  obtain  good  results  at  a patch  that  is 
excited  and  poor  results  at  a patch  that  is  not  excited. 

Another  conclusion  worth  noting  from  all  the  data  presented 
in  figures  3 through  10  is  that  the  plane  wave  illumination 
results  agree  reasonably  well  with  the  code  verification  data 
for  ka  as  large  as  2.3.  In  many  cases  the  agreement  is  still 
reasonable  for  ka=2.9.  These  results  indicate  that  it  is 
possible  to  give  up  a certain  measure  of  accuracy  and  have 
fewer  zones  per  wavelength  than  previously  thought.  For  the 
data  presented,  the  ratio  of  the  wavelength  i to  a zone  dimen- 
sion D is  given  by  A/D=8/ka  which  is  3.5  if  we  accepted  results 
only  up  to  ka=2.3  and  is  2.8  if  we  accept  the  results  up  to 
ka=2.9.  In  either  case  we  see  that  it  is  possible  to  obtain 
acceptable  results  with  fewer  zones  than  has  in  general  been 
previously  thought.  This  can  impact  a scheme  for  determining 
a configuration  of  local  sources.  The  fact  that  sparce  illu- 
mination gave  good  results  also  provides  a rationale  for  employ- 
ing fewer  sources.  Both  of  these  results  can  assist  the  choice 
of  a configuration  to  be  employed  in  an  experiment. 

At  this  point,  it  should  be  noted  that  no  part  of  our  ex- 
plicit sphere  calculation  can  be  used  to  infer  any  experimental 
information  for  very  early  times  since  our  calculation  was  not 
appropriate  for  high  frequencies.  Another  limitation  of  our 
calculation  should  be  pointed  out.  The  sphere  does  not  have 
a sharp  resonance  and  this  could  contribute  to  the  fact  that 
the  patch  containing  the  POE  required  direct  excitation  in 
order  for  good  results  to  be  obtained.  For  structures  having 
more  pronounced  resonances,  it  is  possible  that  near  resonance 
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a given  POE  can  be  excited  without  having  the  source  in  as 
close  proximity  as  indicated  by  our  sphere  results.  Having 
elaborated  on  the  limitations  of  our  sphere  calculations  we 
would  like  to  emphasize  that  the  general  theory  presented  in 
this  work  is  valid  for  all  frequencies  and  consequently  all 
time . 

We  now  address  the  essential  aspects  required  by  our  ana- 
lysis for  each  local  source.  In  the  absence  of  any  other  ob- 
jects and  sources,  their  radiated  fields  should  rapidly  decay 
away  from  their  patch  location  and  at  the  same  time  their 
radiated  fields  should  vary  slowly  over  their  own  patch  loca- 
tion. The  simplest  source  that  possesses  locality  to  some 
extent,  is  a half-loop  placed  above  the  patch.  A simple 
calculation  shows  that  the  fields  decay  rapidly  for  distances 
larger  than  the  radius  of  the  loop.  Despite  its  local  character 
which,  to  a certain  degree,  satisfies  one  of  our  conditions  for 
an  allowable  source, there  are  difficulties  with  the  half-loop 
that  we  will  briefly  discuss:  i)  the  field  due  to  a half-loop 
is  slowly  varying  over  a region  surrounding  the  center  of  the 
loop  but  the  maximum  linear  dimension  of  this  region  is  signi- 
ficantly smaller  than  the  radius  of  the  loop.  To  remedy  this 
we  may  either  consider  a half-loop  much  larger  than  the  patch 
or  a "solenoid"  consisting  of  many  parallel  half-loops  with 
its  dimensions  not  significantly  larger  than  the  dimensions  of 
the  patch.  In  the  case  of  a large  half-loop  the  incident  field 
will  now  vary  rapidly  over  other  patches,  and  we  have  not 
assessed  the  effect  of  this  behavior  in  our  calculations. 
However,  the  numerical  solution  is  only  a convenience  for 
studying  selective  patch  excitation  and  its  inapplicability 
does  not  invalidate  the  potential  use  of  the  half-loop  as  a 
portable  simulator.  ii)  The  "solenoid"  is  an  improvement  with 
regard  to  the  condition  of  slow  variation  but  it,  as  well  as 
the  large  half-loop,  may  interact  with  the  sphere  substantially 
and  this  could  significantly  alter  the  transfer  operator 
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as  explained  earlier  in  connection  to  non-rigidity  of  sources. 
Despite  all  the  described  limitations,  both  the  "solenoid" 
and  the  half-loop  have  sufficiently  desirable  features  to  be 
included  in  an  experimental  program. 

Finally,  we  discuss  the  Singularity  Expansion  Method  (SEfl) 
as  it  relates  to  alternate  simulation.  We  do  this  because  it 
offers  a hope  of  determining  the  global  capabilities  of  a 
configuration  of  portable  sources.  We  will  now  interpret  our 
results  as  related  to  SEM.  An  SEM  external  interaction  solution 
has  the  form 


J(r,y)  \(v)  ^ 

The  natural  modes  v and  natural  frequencies  > c are  intrinsic 

—a  a 

properties  of  the  metallic  body.  The  coupling  coefficients 
na  depend  on  both  the  coupling  vectors  (also  an  intrinsic 
property  of  the  body)  and  the  incident  field.  Thus,  once  the 
natural  modes,  the  coupling  vectors,  and  the  natural  frequencies 
are  known,  the  responses  to  various  excitations  in  the  SEM 
prescription  are  obtained  by  determining  the  corresponding 
coupling  coefficients. 

Admittedly  there  is  no  known  recipe  for  obtaining  the 
coupling  coefficients,  in  general,  but  at  least  we  know  that 
for  the  sphere  and  plane  wave  illumination  the  correct  coupling 
coefficients  are  class  1 given  by 


"nn'mop’''1  = e 


^nn,-1r)cto 


J , J 
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me 
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dvH. 
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where 
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Y , are  the  pole  locations  (y  ,c  = natural  frequencies),  t~ 

is  the  instant  at  which  the  incident  wavefront  hits  the  sphere, 
i nc  ^ i. nc 

J =nxH  and  p stands  for  polarization,  A are  the  eiqen- 
— p -p  1 n J 

values  of  the  Magnetic  Field  Integral  Operator  L,  and  J are 

^ r nmo 

eigenfunctions  of  L corresponding  to  eigenvalue  1-A^. 

If  we  were  to  compare  responses  to  selective  patch  excita- 
tion and  plane  wave  illumination,  we  could  assume  that  the 
coupling  coefficients  for  patch  excitation  are  also  given  by 
equation  1 and  proceed  to  calculate  them.  The  comparison  of 
the  coupling  coefficients  for  the  two  excitations  would  allow 
us  then  to  ascertain  how  well  selective  patch  excitation 
simulates  plane  wave  illumination.  At  this  point,  however, 
caution  should  be  exercised.  To  clarify  the  point  we  are 
trying  to  emphasize,  consider  the  case  whereby  we  excite  all 
patches  on  the  sphere  but  one,  in  the  manner  that  was  explained. 
The  MFIE  solution  shows  that  the  total  current  induced  on  the 
sphere  is  everywhere  approximately  equal  to  the  current  for 
plane  wave  illumination  except  at  the  center  of  the  patch 
that  was  not  excited.  However,  if  we  were  to  use  SEM  for  the 
comparison  of  the  two  types  of  excitation,  the  coupling 
coefficients  for  the  first  few  modes  would  be  approximately 
equal  and  this  result  might  lead  one  to  the  false  conclusion 
that  the  simulation  was  adequate.  Notice,  however,  that  our 
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patch  zoning  results  provide  no  information  as  to  any  early- 
time  SEM  results  and/or  conclusions. 
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